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ion collisions are shortly reviewed. The common features of the models, responsible for 
^ ■ the successful description of hadronic observables, are outlined. 

in 

I. INTRODUCTION 

In this talk I present several models which turned out to be quite successful in repro- 
ducing the experimentally measured hadron spectra in nucleus-nucleus collisions. The 
discussed models use concepts borrowed from relativistic hydrodynamics but they do not 
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Several popular parameterizations of the freeze-out conditions in ultra-relativistic heavy- 



include the complete time evolution of the system. For this reason, they may be called 
hydro-inspired (freeze-out) models. The measured particle spectra reflect properties of 
matter at the stage when particles stop to interact. This moment is called the kinetic 
(thermal) freeze-out. Hydro-inspired models help us to verify the idea that matter, just 
before the kinetic freeze-out, is locally thermalized and exhibits collective behavior, such 
as the transverse and longitudinal expansion. If this is really the case, i.e., if the hydro- 
inspired models describe the data well, we may infer the thermodynamic properties of 
matter at freeze-out, such as the values of the temperature and flow, and request that the 
advanced hydrodynamic models reproduce this configuration. 

In my opinion, the real aim of the freeze-out models is to form a simple link between 
sophisticated hydrodynamic calculations [ HJ 121 El IH El E] describing the full time evo- 
lution of matter (with the inclusion of the phase transition) and the rich bulk of the 
experimental data describing soft phenomena. This is an appealing idea, however, one 
encounters several problems on the way to achieve this task, since certain features of the 
hydrodynamic models and the freeze-out models are quite different. For example, at the 
first sight one can realize that typical shapes of the freeze-out hypersurfaces used in the 
advanced hydro calculations and in the freeze-out models are quite different. Clearly, 
further work is necessary to make the connection between the hydro calculations and the 
freeze-out models more consistent. 



*This work was supported in part by the Polish State Committee of Scientific Research, Grant No. 2 
P03B 05925. 
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In any case, an attractive feature of the hydro-inspired models is that they are very 
effective parameterizations of the final state, that use few parameters possessing clear 
physical interpretation. As long as we do not all have hydrodynamic codes that we are able 
to run on our PCs, the freeze-out models form a very convenient and easy accessible tool 
to interpret the data. In this talk I intend to review several models which give consistent 
description of many observables (not only of the particle spectra). Their list includes: 
different versions of the blast- wave model [HE], the Buda-Lund model [ H3 EH E] , the 
Seattle model [11211121, the Durham model [CUEH], the Cracow single-freeze-out model 
[ EH HH EU , and THERMINATOR [ H2| ■ 

Below, I will refer only to the Au+Au collisions studied at the top RHIC energies; I 
want to show how different models describe the same set of data rather than to show how 
one model is able to describe different physical situations. Examples of the application of 
the thermal approach to calculate the spectra at lower energies may be found in Refs. [ 

mino- 

2. COOPER-FRYE FORMULA AND EMISSION FUNCTION 

Before going into the discussion of the blast-wave model, I would like to make a few 
remarks about the Cooper- Frye formula [ I22| . It is frequently used in the hydrodynamic 
calculations to obtain the momentum distribution of the emitted particles. In this talk, 
the Cooper- Frye formula will serve us as a reference point for the characteristics of different 
freeze-out models. Its standard form is 

EP % = l P^to/eq (?■«(*))> C 1 ) 

where E p = y/ m 2 + p 2 is the energy of a particle, c?S M (x) a three-dimensional element of 
the freeze-out hypersurface, u^(x) is the hydrodynamic flow, and / eq is the equilibrium 
distribution function. In the general case, the hypersurface includes the time-like (defined 
here by the condition: cE^dE^ = (c/S ) 2 - (d£i) 2 - (rfS 2 ) 2 - (d£ 3 ) 2 > 0) and space- 
like (<iX M <iX M < 0) parts. The advanced hydrodynamic calculations include both parts, 
while hydro-inspired models include typically only the time-like parts. The success of the 
hydro-inspired models may indicate that the contributions from the space-like parts are 
negligible. On the other hand, detailed microscopic studies show that the contributions 
from those parts are substantial and difficult to include in the consistent way, unless one 
uses the approach based on the kinetic theory [ EHl ESI EEl EZl EHI EHl E01 EH E21 ESJ - 

The Cooper-Frye formula may be rewritten in such a way that the distribution of the 
particles in the momentum space is given as the space-time integral over the so-called 
emission or source function 

= l dAX l P^Mx'W - X ) ^ (P • <X')) EE j d 4 xS(x, P ). (2) 

The emission function defines the space-time distribution of the points from which the 
observed hadrons are emitted. Very often Eq. (J2J) is used with S(x,p) modeled without 
any reference to the Cooper-Frye formula. Such a procedure lays emphasis on the fact 
that the hydro-inspired models should aim mainly at the reconstruction of the realistic 
emission function. 
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typical freeze-out curve in 
hydro-inspired (freeze-out) 
models 



z = 



typical freeze-out curve in 
hydrodynamic models 



Figure 1. Freeze-out curves in the Minkowski space, i.e., projections of the freeze-out 
hypersurfaces on the plane z = for cylindrically symmetric systems, in the advanced 
hydro calculations (dashed line) and in the freeze-out models (solid line). 



3. BLAST- WAVE MODELS 

Different versions of the blast-wave model originate from the paper by Siemens and 
Rasmussen [EH, where a relativistic formula for the particle distribution corresponding 
to a thermalized and radially expanding system was first given. More recent applications 
use the same concepts but different geometries of the expansion, more suitable for the 
description of the ultra-relativistic heavy- ion collisions, are considered [ EE] . 

3.1. Cylindrically symmetric systems 

For boost-invariant and cylindrically symmetric systems, the Cooper-Frye formula 
leads to the very popular model of Schnedermann, Sollfrank and Heinz [|7j. For constant 
transverse flow, vt = tanhp = const, one obtains the rapidity and transverse- momentum 
distribution of the emitted particles in the form 

dN } dr 

-m T Ki [/?m T cosh(p)] I Q [Pp T smh(p)} / d( r(()t(()- 



dyd 2 pT 2ir 2 J d( 

~^PtK [/5m T cosh(p)] h [(3p T smh(p)] J dQ KO*(C)^, (3) 



where wit = y + p\ is the transverse mass, (3 = 1/T is the inverse temperature, p is the 



chemical potential, p is the transverse rapidity, t and r = ^Jr 2 + r 2 are the coordinates of 
the freeze-out hypersurface at z = 0, while _K"o,i and Jo,i are the modified Bessel functions. 
The integrals on the right-hand-side of Eq. Q involve the parameterizations of the freeze- 
out position coordinates and times in terms of the parameter (. In practical applications, 
the second line on the right-hand-side of Eq. (jHJ) is very often neglected, which means 
that the variations of the emission times are ignored. This procedure implicitly denotes 
that the particle emission takes place at a constant laboratory time (at z = 0), see Fig.[TJ 
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An example of the procedure outlined above is the blast-wave fit performed recently 
by the BRAHMS Collaboration [EH]. The optimal value of the temperature obtained 
from the fit to the Au+Au data (^/s^n = 200 GeV, midrapidity region, centrality class 
0-10%) is 110 MeV, while the average transverse flow velocity equals 0.65 c. In this, and 
other similar cases, we have to be aware that a theoretical boost-invariant model has 
been applied to the system which is not boost-invariant as a whole. In my opinion, such 
procedure is justified only at the top RHIC energies in the midrapidity region where the 
rapidity distribution may be regarded as flat within 1 or at most 2 units of rapidity. 

3.2. Cylindrically non-symmetric systems 

For boost-invariant and cylindrically non-symmetric systems one may use the param- 
eterization of the emission function introduced by Retiere and Lisa It takes into 
account the possible ellipsoidal shape of the system created in non central collisions. The 
emission function proposed in [[Sj may be obtained from the Cooper- Frye formula if the 
constant evolution time, r = tq = const, is replaced by a gaussian distribution 

~(^-tq) 2 1 

S (x, p) = Zm T cosh(r? - y) ft(r, <f> a ) e ^ T - , (4) 

Here Z is an arbitrary normalization constant, r\ — 1/2 In (t + z)/(t — z) is the spacetime 
rapidity, tan0 s = r y /r x , At is the emission time, and Q describes the spatial distribution 
of matter in the transverse plane 



Q(r, <f> a ) = Q(r) = - - , r(r, <f> a ) = 

The parameter a s describes the surface diffuseness of the source. In the limit a s — ► 0, the 
matter in the transverse plane is confined in an ellipse defined by the parameters R x and 
R y . The explicit form of the flow field (and of the expression p ■ u) is given in Ref. [ 
|H] , here we only note that the expansion velocity in the transverse plane is perpendicular 
to the elliptical shell confining the matter. 

The optimal values of the temperature obtained by Retiere and Lisa in their analysis 
of three different centrality classes (0-5%, 15-30%, 60-92%, the Au+Au collisions at the 
beam energy of = 130 GeV, combined PHENIX and STAR data [EH EH]) are: 

107, 106, and 100 MeV, respectively. The corresponding values of the transverse flow are: 
0.52 c, 0.50 c, and 0.47 c. 

3.3. Summary on the blast-wave models 

When applied to the Au+Au collisions at the top RHIC energies, the blast-wave models 
give good description of the p^-spectra, the elliptic flow, and the HBT radii at midrapidity. 
An interesting feature of the Retiere-Lisa model is that it includes, as the special case, 
the situation where the in-plane flow is not stronger than the out-of-plane flow, but the 
correct sign and magnitude of the coefficient v-i is reproduced. This feature is realized 
by the geometry of the model which implies in this case that more matter flows in the 
reaction plane than out of the reaction plane. 

The blast- wave fits to the data indicate short evolution times, r < 10 fm, and very 
short emission times, At < 1 fm. The typical value of the temperature is about 100 



(rcos(0 s )) 2 (rsin(0 s )) 2 



(5) 
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MeV. The blast-wave models do not predict the absolute normalization of the spectra, 
hence an extra parameter is required to normalize each considered spectrum. In addition, 
many of the applications of the blast-wave model do not include the effects of the decays 
of resonances. 



4. BUDA-LUND MODEL 

The standard version of the Buda-Lund model was formulated to describe cylindrically 
symmetric systems with no constraint of the boost-invariance [ EH . Later the model 
was extended to describe ellipsoidally symmetric systems [ ■ The standard emission 
function of the model has the form 



q( \ _9 m T cosier? - y) 1 

(2tt)3 exp _ ± ! (27rAr»)V» ex P 



2Ar 2 



(6) 



where the temperature and chemical potential depend on the position coordinates 

T{x) " 
1 



T(x) 




(7) 



(8) 



In Eq. (jHJ), the quantity T is the temperature in the center at the mean freeze-out time, T s 
is the temperature on the surface at the mean freeze-out time, and T e is the temperature 
in the center at the end of particle emission. The flow pattern assumed in the Buda-Lund 
model has a Hubble-like structure, with the velocity of the fluid element proportional to 
its distance from the center. Such patterns appear in the analytic and numerical [ |4*U] 
solutions of the equations of the relativistic hydrodynamics. For physical interpretation 
of other model parameters we refer to the original papers [19*1 I10|. 

The Buda-Lund model gives a very good description of the pr-spectra, t> 2 and HBT 
radii in the full rapidity range. The fits to the experimental data indicate small values 
of the evolution and emission times which are compatible with those obtained from the 
blast-wave analysis (for Au+Au at ^snn = 200 GeV the optimal Buda-Lund fit yields: 
T = 196 MeV, T e = 117 MeV, T s = 89.7 MeV, R G = 13.5 fm, R s = 12.4 fm, r = 
5.8 fm, At = 0.9 fm, and At] = 3.1). The effects of the decays of resonances are not 
included explicitly and the information about the normalization of the spectra is derived 
from the HBT data. The Buda-Lund fits indicate a freeze-out temperature distribution, 
with an average freeze-out temperature of 100 - 120 MeV. However, about l/8th of the 
particles is found to be emitted from a very hot center with T ~ 200 MeV. A comparison 
of this result with the lattice QCD results for the critical temperature was considered as 
an indication of quark deconfinement in Ref. 



5. SEATTLE MODEL 



The Seattle model formulated in Refs. [ EI] describes cylindrically symmetric sys- 
tems. No assumption about the boost-invariance is made. The main characteristic feature 
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of the model is inclusion of the effects of the final-state interaction of the outgoing pions. 
Such effects are taken into account by the use of the distorted wave functions ty^'(x). 
The formalism, following Ref. [ I42| . is based on the generalized emission function, which 
may be used to get the one-particle and two-particle distributions in the momentum space 

S{x,p,q) = Jd 4 K'S (x,K') l0^ e - iK '-^(x + xy2)¥ p2 >(x-x'/2), (9) 

p = (p 1 +p 2 )/2, q = Pi-P2- (10) 

The quantities p\ and P2 are the pion four-momenta. In the special case, where the 
distorted wave functions are replaced by the plane waves, the standard formulation is 
recovered with the emission function reduced to So, 

¥-\x)^e^ x , S(x,p,q)^S (x,p)e i ^, S(x,p, 0) - S Q (x,p). (11) 

The function So resembles the parameterization used in the Buda-Lund model, however, 
in this case the thermodynamic parameters are constant, 



m T cosh rj 
S (x,p) = . exp 
/2tt(At) 2 



[t - r ) 2 r? 2 



2Ar 2 2A^ 2 



n(r) 1 , . 

— — . (12) 

(27r) 3 exp(^ f ^) - 1 1 ' 



The distribution of matter is characterized by the function 
1 



n(0 = 7 — 7—7-^ — v r = v r ' + i- ( 13 ) 

hp ) + !] 

The optimal values of the parameters obtained from the fit to the STAR data [ 03] 
describing the Au+Au collisions at y77\7v = 200 GeV are the following [EH]: T = 215 
MeV, fa = 123 MeV, r = 8 fm, At = 2.7 fm, R ws = 12 fm, a ws = 0.8 fm, Ar] = 
1, 77/ = 1.5, w = 0.142±0.046 fm" 2 , and w 2 = 0.582 ± 0.014 + i(0.123±0.002). The 
parameter rjf defines the magnitude of the transverse flow [H2]) whereas the parameters 
w and w 2 define the optical potential which modifies the outgoing pion wave functions. 
The strength of the attraction inside the medium is greater than m 2 ,, hence a pion behaves 
to large extent as a massless particle. This behavior may be related to the phenomena 
discussed in Refs. [ EH E3] . 

The Seattle model describes successfully the pion transverse-momentum spectra and 
the HBT radii in the region px > 100 MeV. On the other hand, the model predicts non- 
monotonic structures in the transverse-momentum region below 70 MeV. Such predictions 
were confronted with the PHOBOS data in Ref. [ |4"B] and certain discrepancies between 
the data and the model were found. Very likely the non-monotonic structures are due 
to the sudden change of the optical potential in the transverse direction and a better 
agreement with the data may be achieved by the modification of this dependence. 

6. DURHAM MODEL 

The Durham model, proposed by Renk in Refs. [El E3; is based on the parame- 
terization of the hydrodynamic expansion in the proper-time interval tq < t < tj. The 
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pressure profiles determine the evolution of the boundaries of the system, then the bound- 
aries determine the evolution of the entropy density. In the next step, the equation of 
state is used to calculate the pressure. The iterative procedure ensures that this approach 
is self-consistent, i.e., the pressure profiles which determine the evolution of the system 
agree with the pressure obtained from the equation of state. The freeze-out happens at a 
fixed temperature Tf and the Cooper-Frye formula is used to calculate the spectra. The 
optimal parameters obtained from the fit to the PHENIX data [ ^\ HH] include: r = 
1 fm, Tf = 19 fm, Tf — 110 MeV. The Durham model describes nicely the p^-spectra and 
HBT radii. No resonance decays are included in the calculation of the spectra. In my 
opinion, a direct comparison with the real hydrodynamic calculations will strengthen the 
importance of this model, if the validity of the proposed parameterization is confirmed. 



7. CRACOW SINGLE-FREEZE-OUT MODEL 

The standard version of this model [ El describes boost-invariant and cylindrically 
symmetric systems. The model is based on the extreme assumption that the chemical 
and kinetic freeze-outs coincide. The thermodynamic parameters of the model are ob- 
tained from the analysis of the ratios of hadron abundances [ HH] . In this respect the 
Cracow model is closely related to different statistical models developed recently by many 
groups [Eni EU E2 E3 EH ES]- An important feature of the model is that a complete 
set of hadronic resonances is included in the calculation of various physical observables. 
The resonances are included by the use of the emission function that corresponds to 
the Cooper-Frye formula convoluted with the momentum splitting functions B and the 
space-time displacement functions 5, 

S {x 1 ,p 1 ) = E 



d 3 pid i Xi 

d?> V2 „ _ v f ,_ „ „-r,r, f j4_ cU) (_ , P2T2 



/ ^B(p 2 , Pl ) J dr 2 T 2 e- r ^J d A x 2 5^ [x 2 



J P2 

X 



m 2 



dT,^ (x N ) p% <5 {4) ( x N + PnTn - x N A f N [p N ■ u (x N )} (14) 
V m N J 



Here the indices 1,...,N denote particles in a cascade of decaying resonances (N refers 
to the resonance decaying on the freeze-out hypersurface, while 1 refers to the final ob- 
served particle, e.g., a pion, other symbols are defined in H3 EE])- The freeze-out 
hypersurface is defined by the conditions: 



t 2 — r 2 x — r 2 — r\ = const, r = ^r 2 + r 2 < r max (15) 
u ^ = ^ = _L (l^JjL^\ dW = u» Tf nv sinh(p)cosh(p) dpdndcf). (16) 

7"inv ^inv V tut/ 



The standard version of the model uses only 4 parameters: temperature (T), baryon chem- 
ical potential (/xb), T inv , and r max . The values of the temperature and baryon chemical 
potential, obtained from the analysis of the ratios of hadron abundances, are independent 
of centrality. For the top RHIC energies one finds T = 165.5 MeV and p# = 28.5 MeV. 
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The expansion parameters r inv and r max are determined by the fits to the pr-spectra, and 
their values depend on centrality 

The Cracow model gives a good description of the ratios of hadronic observables, pt- 
spectra, the HBT radii -R S id e and -R ou t- It also gives satisfactory description of the balance 
functions [ EZI and of the invariant masses of the pion pairs [ EH] • Generalizations of 
the model describing systems without cylindrical symmetry [EH] or boost-invariance [ l3*o] 
were proposed to account for the observed values of the coefficient v 2 and the pr-spectra 
at finite values of the rapidity as measured by BRAHMS [ 60 . 

8. THERMINATOR, SHARE, THERMUS 

THERMINATOR [ 19j is the Monte-Carlo version of the Cracow and blast-wave models 
(the blast-wave model has been extended in this case to include the complete set of 
hadronic resonances). The program uses the same input from the Particle Data Tables 
[ IBT] as SHARE [E2|. With the physical input of the thermodynamic and expansion 
parameters, THERMINATOR delivers the full space-time information about cascades of 
decaying resonances. The values of the thermodynamic parameters may be taken from 
SHARE or from other programs used to study the ratios of hadronic abundances, for 
example from THERMUS [63J. The information about the space-time positions of the 
produced hadrons may be used to study different types of correlations. Moreover, the 
Monte-Carlo approach facilitates the inclusion of the experimental cuts and acceptance. 

9. HUMANIC'S RESCATTERING MODEL 

In the end of this talk I would like to note that simple transport models may play a 
similar role as hydro-inspired models. An example of such approach is Humanic's model 
[ lb'4| . This model is a Monte-Carlo simulation of the evolution of purely hadronic system. 
The initialization stage may be considered as equivalent to the hadronization process. The 
transverse geometry is determined by the overlapping region of the two colliding nuclei, 
while the initial p?- and ^-distributions are assumed to have the following shape 

1 dN _ Q mT dN 

rriT drriT exp(rriT/T) ± 1 ' dy 

with the optimal parameters: T = 300 MeV and a y = 2.4. The longitudinal positions and 
times of the initialized hadrons are: z^&d = Thad sinhy and thad = Thad coshy, where Th a d is 
the initial (proper) time, with the standard value of 1 fm. 

The subsequent evolution of the system includes binary collisions of: tt, K, N, A, A, 
p, u>, T], rj', (ft and K*. The freeze-out happens (dynamically) at about 30 fm (Au+Au, 
130 GeV). The model gives successful description of the slope parameters, V2(m,pT,rj), 
and the p? and centrality dependence of the HBT radii. The results are sensitive to t^^] 
larger values of r had imply fewer collisions and the rescattering-generated flow is reduced. 

10. CONCLUSIONS 

A common feature of the hydro-inspired models is the short evolution time, r < 10 fm, 
and even shorter emission time At << t . Another common feature is a large value of 
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the transverse flow, (v t ) ~ 0.5 c. Altogether, these observations indicate an explosive 
scenario for the Au+Au collisions observed at the top RHIC energies. 

Models which do not include resonances or pion final-state interactions, yield rather low 
values of the freeze-out temperature, T ~ 100 MeV. We know, however, that the effects 
of the resonance decays must be important, since large abundances of certain resonances 
were observed [ EE] • In this situation, I think that the hydro-inspired models should 
enter a new stage where the effects of both the transverse flow and the resonance decays 
are commonly taken into account. Another development of the hydro- inspired models 
should include the emission from the space-like parts. An example of such approach is a 
recent paper by the Kiev-Nantes group [KT7j. 
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